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Introduction 

Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in 
hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other 
blades. This interaction, known as blade-vortex interactions(BVI), results in unsteady loading ot the blades 
and can cause a distinctive acoustic signature. 

Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interac- 
tions can help rotor designers and engineers to predict rotor performance and to develop designs for low 
acoustic signature. Such a predictive method must preserve a blade’s shed vortex for several blade revolu- 
tions before being dissipated. A number of researchers have explored the requirements for tills task. 

Tliis paper will outline some new capabilities that have been added to the NASA Ames’ OVERFLOW 
code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these 
improvements, a number of case studies will be presented. These case studies consist of free convection ot a 
2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a lull 3-D unsteady viscous 
solution of a helicopter rotor in forward flight 

In tliis study both central and upwind difference schemes are modified to be more accurate. Central dif- 
ference scheme is chosen for tliis simulation because the flowfield is not dominated by strong shocks. The 
feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interac- 
tion. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in tully- 
implicit manner at each time-step. The spatial accuracy is cither second and fourth-order central difference 
or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and dem- 
onstrates the methods for several sample cases and for a helicopter rotor. Preliminary computations on a 
rotor were performed by using tliis method! 1 ] and is in the process of documentation. 

Model Problems 

The scheme is tested by calculating several model problems involving vortex convection, moving, and 
oscillating grids. The method also provides the various aspect of the moving grid computations. Two cases 
presented for this study are closely related to typical helicopter flowfields. Recommendations will be pre- 
sented for practical implementation and use of the various numerical schemes. 

The first test case is a passive convection of vortex in a free stream and monitoring a measure ot the rate 
of decay of the vortex. A perfect method should convect the vortex in space with no dissipation or disper- 
sion. This case tests the time-accurate free convection of a vortex flow. Consider a compressible Euler equa- 
tions of gas dynamics problem for a 2D rectangular domain. The mean flow is p - 1, p- 1 , and 
(t(,v)= (0,0) (steady flow)or («,v)- (1,0) (horizontal flow), or (», v)= (1,1) (diagonal flow). We add, to 
the mean flow, an isentropic vortex i.e., perturbation in («, v) and the temperature T = p/p , no perturbation 
in entropy. The computational domain is taken as [0, 10 ] x [0, 101, extended periodically in both directions. 
The convection speed is Mach 0.5. The exact solution of the Euler equation with tltese initial and boundary 
conditions is simply the passive convection of the vortex with the mean velocity. Periodic boundary condi- 
tion allows us to compute the long time evolution ot die vortex. No shocks arc considered. 

The initial vortex is positioned at tlx; center of the domain. The grid resolution is 8 1 x 8 1 and the 
timestep DT is normalized by the freestream speed of sound and the imit of length ol the domain. II DT = 


1 


o.l , it takes an elapsed time of 20 for the vortex to convert through one period and (ideally) return back to 
the center of the grid. 

Various methods were tested with this problem. Plots arc shown tor the density variable. The initial vor- 
tex center is represented with grid lines. Tliis facilitates showing the phase shill ol the vortex alter time evo- 
lution. Fig. I shows the contour plots for the first period of the vortex traverse. Initial start-up vortex is shown 
in Fig 1 .(a). Fig. l(b-c) is for first-order time accurate computation with and without subiteration respec- 
tively. Fig. 1(d) is with the second order in time and wilh three subiterations, and for the first period, this 
closely resemble the initial vortex. For a bigger time-step, and nosubiteration (Fig lc) the vortex is almost 
out of phase at the very first period. The rest of the line plots of Fig. 2 thorough 1 2 show the density at a cut 
section of z = 5 through the domain for the various methods and options used. 

In the second problem, a case study of unsteady 2D flows over an oscillating NACA0015 airfoil have 
been performed for an attached flow. The results were compared with experimental data as well as with the 
existing various numerical computations^]. The experiment chosen for all tlie previous computations were 
done at NASA Ames Research Center. The experiment gave detail measurements of static and dynamic 
pressure distributions and cycle averaged lilt, drag, and pitching moment coefficients. 

Fig. 13 shows hysteresis of lift, drag, and pitcliing moment coefficients for attached flow. The only tur- 
bulence model used in tliis study is Baldwin-Lomax model. Computations were done on a C-type grid of 
total of 259 x 60 points, of which 100 points are on the upper surface, 80 points on the lower surface of tlx: 
airfoil, 40 points along the wake cut. The farfield boundary extended to 15 chord length away from the air- 
foil surface. In the present study, only 1440 steps were taken for each pitcliing cycle, yet comparable accu- 
racy were obtained with the central difference scheme of both second and fourth-order. Third-order upwind 
MUSCLE scheme shows considerable discrepancies in the drag and pitching moment prediction. 

Results of third-order scheme of Fig. 1 3 were obtained without any limiter to the MUSCLE scheme. Fig. 
14 compares third-order MUSCLE scheme results with and without limiters. All of these calculations agree 
well with the computations by other investigators. From these calculations, it is evident that no limiter 
should be used for flow without shocks and for cases with shocks, choice of proper limiter is very important. 

Helicopter Rotor Problem 

This study will demonstrate a method of computing flowfield of arbitrary rotors in forward flight. The 
equations are solved on a system of overset grids that allow for prescribed cyclic and flapping blade motions. 
Computed results will be compared either with flight test data or wind-tunnel experiment with a scaled rotor. 
A full 3-D simulation with rotors showing promising results were shown in Ref. 1. Details of this will be 
presented in the full paper. 


Bibliography 

1. Strawn, R.C., Ahmad, J., Duque, E.P.N., “Rotorcraft Aeroacoustics Computations with Overset-Grid CFD 
Methods,” 54th Annual Fontm, American Helicopter Society, May 22-24, 1998, Washington, D.C. 

2. Sungho Ko and W. J. McCroskey, “Computation of Unsteady Separating Flows over an Oscillating Air- 
foil,” AIAA-95-0312, 33rd Aerospace Sciences Meeting and Exhibit, January 1995, Reno, NV. 


2 



3 






Figure 4. Effect of number of subiterations. First order time accurate. 
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Figure 5. 1st order time accurate. Central difference 4th and 2nd order spatial accuracy with tourth order dissipation and third order MUS- 
CLE scheme with ROE upwinding is with limiter turned on. 
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Figure 6. 2nd order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order 
MUSCLE scheme with ROE upwinding is with limiter turned on. 
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Figure 7. 1st order time accurate. Central difference 4th and 2nd order spatial accuracy with tourth order dissipation and third order MUS- 
CLE scheme with ROE upwinding is with limiter turned off. 
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Figure 8. 2nd order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order 
MUSCLE scheme with ROE upwinding is with limiter turned off. 
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Figure 9. 1st order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order MUS 
CLE scheme with ROE upwinding is with limiter turned on. 
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Figure 10. 2nd order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order 
MUSCLE scheme with ROE upwinding is with limiter turned on. 
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Figure 1 L 1st order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order 
MUSCLE scheme with ROE upwinding is with limiter turned off. 
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Figure 12. 2nd order time accurate. Central difference 4th and 2nd order spatial accuracy with fourth order dissipation and third order 
MUSCLE scheme with ROE upwinding is with limiter turned off. 








